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A model of discrete dynamics of entanglement of bipartite 
quantum state is considered. It involves a global unitary dy- 
namics of the system and periodic actions of local bistochastic 
or decaying channel. For initially pure states the decay of en- 
tanglement is accompanied with an increase of von Neumann 
entropy of the system. We observe and discuss revivals of 
entanglement due to unitary interaction of subsystems. For 
some mixed states having different marginal entropies of the 
subsystems we find an asymmetry in speed of entanglement 
decay. The entanglement of these states decreases faster, if 
the depolarizing channel acts on the "classical" subsystem, 
characterized by smaller marginal entropy. 

PACS numbers: 03.65 Bz, 03.67.-a 



ment. We demonstrate a possible presence of revivals of 
entanglement caused by the global unitary evolution en- 
tangling the subsystems between consecutive actions of 
the environment. 

The paper is organised as follows. In section II we 
describe a simple model of discrete time evolution. In 
section III we derive bounds on the entropy increase un- 
der the action of the environment. Then in section IV 
we analyse the decrease of entanglement versus increase 
of the degree of mixing of the initially pure states. The 
asymmetry in the entanglement decay depending on the 
subsystem subjecting to influence of environment is de- 
scribed in section V. In section VI we consider the entan- 
glement revivals. The results of the paper are discussed 
in section VII. 
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I. INTRODUCTION 

Quantum entanglement is one of the most subtle and 
intriguing phenomena in nature Its potential use- 

fulness has been demonstrated in various applications 
like quantum teleportation, quantum cryptography and 
quantum dense coding. On the other hand, quantum en- 
tanglement is a fragile feature, which can be destroyed 
by interaction with the environment. This effect due to 
decoherence ^ , is the main obstacle for practical imple- 
mentation of quantum computing. A model allowing to 
study the dynamics of entanglement in presence of inter- 
action with the environment has been recently analysed 
by Yi and Sun H). 

In this paper we investigate destruction of the entan- 
glement in a proposed model of discrete dynamics. We 
consider a simple bipartite system consisting of two spin- 
^ particles. Only one of them is subjected to periodic 
actions of a quantum channel, which represents the inter- 
action with environment. As the initial states we choose 
random states taken from the ensemble of pure separa- 
ble states and from the ensemble of maximally entangled 
pure states. We also investigate the time evolution of 
mixed states having some special property. The corre- 
sponding system is composed of two subsystems exhibit- 
ing different properties with respect to some entropy in- 
equality which is satisfied by all classical systems. One of 
the subsystems satisfies the inequality and may be con- 
sidered "classical", while the other, "quantum" subsys- 
tem violates the inequality. We investigate an asymme- 
try in the process of destruction of entanglement with 
respect to the subsystem interacting with the environ- 



II. MODELS OF TIME EVOLUTION 

In this paper we consider the bipartite state subjected 
sequential interactions with environment. They are mod- 
elled by quantum channels, defined as completely positive 
linear maps, preserving the trace of the state ||]. 

Let a he a density operator acting on a finite- 
dimensional Hilbert space Ti.. The most general form 
of the quantum channel is the following transformation 
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a' = A(cr) = ^ F.ctF/, where 



i=l 



I (1) 



If in addition ^i^/ — holds then the channel is 

called bistochastic. 

Bistochastic channels can be alternatively defined as 
channels which do not decrease von Neumann entropy of 
quantum states. 

A particular example of the bistochastic channel is 
given by random external fields . They can be writ- 
ten as 
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, K are unitary operators and the 



where Ai, i — 1,2, 
vector of probabilities p = [pi, . . . , px] is normalised 
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^p^ = l, Pi > 0. 
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Such random systems can be described in the formahsm 
of quantum iterated function systems 0. The so called 
Kraus form (|]) can be reproduced setting Vi = ^fplAi. It 
is worth to note that in the case of the most elementary 
quantum system described on the Hilbert space Ti. ~ 
the channel is bistochastic if and only if it is a random 
external field (|^) (see Q ) . Note that an unitary evolution 
of the system can be considered as the simplest case of 
the bistochastic quantum channel with K = 1. 

There exist, however, many quantum channels which 
are not bistochastic. We shall consider the following de- 
caying channel, sometimes called |^ the amplitude damp- 
ing channel) 



a' = Kd{<t) ^ MicrMi + K'hcrNh 



(4) 



where the matrices Mi 
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and A'h — 



are written in the standard basis. 







Let Q denote a mixed state of a 2 (g) 2 system i. 
e. the density operator defined on the Hilbert space 
H = Ha^T-Lb = (E>C'^. The system consists of two 
subsystems A and B which can represent spin-i particles 
or two-level atoms. In our model the unitary dynamics 
is interrupted by periodic actions of the environment as 
shown schematically in Fig. 1. 

Discrete time evolution of the state g reads in our 
model 



g{n+l) = Ug'{n)U^ = U{A{g{n))U'' 



(5) 



where A — I ^ A and the channel A is either bistochas- 
tic (H) or decaying (^. Here U = e*"^ represents a 
unitary transformation which involves an interaction be- 
tween both subsystems A and B described by the Hamil- 
tonian H. We use the dimensionless units and a stands 
for a coupling parameter. Subsequently we shall con- 
sider the cases with H equal either to a j.® ay = H or to 



H' 



In general we shall use four types of dynamics defined 
by four different operators A-s in the formula (||). Three 
of them will be random external fields A^ (||) , all defined 
by the same set oi K = A unitary operators: Ai — I, 
^2 = (Ti, A3 = (72, ^4 = CTs (where ai denote Pauli 
matrices), but with different vectors of probability: (^: 
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e,0,0,e] 
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e, -, -, - , < e < 1. 
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Each dynamics depends on two continuous parameters: 



a contained in U 



governing the unitary dynamics 



and e, included in the vector of probabilities, and describ- 
ing the strength of the coupling with the environment. 
Additional discrete index j labels the different vectors of 



probability, . For these three models of dynamics we 
shall use the compact notation ^. The fourth dynam- 
ics denoted by Qa.p is defined by putting in formula (||) 
the decaying channel (^. Dynamics involving the op- 
eration U with "reflected" (i. e. obtained from H by 
permutating subsystems) Hamiltonian H' = ay® ax will 
be denoted by the same symbols with only one change: 

e. 

Remark.- If a is equal to zero, then the unitary oper- 
ation U in (|^) is reduced to identity transformation. In 
particular, it can be seen that the dynamics Oq ^ corre- 
sponds to periodic action of depolarizing channel . 

Now the essence of our study is the following: we con- 
sider composite quantum systems subjected to the local 
interaction with the environment, which acts on one sub- 
system only. We investigate, how the decay of the entan- 
glement in the system depends on the initial state and 
the type of the dynamics. In particular we analyse, to 
which extent the decrease of the mean entanglement is 
reflected by the evolution of von Neumann entropy of the 
system. 



III. BOUNDS ON ENTROPY INCREASE UNDER 
LOCAL CHANNEL 

We start establishing bounds for the increase of von 
Neumann entropy. 

Proposition .- Under a local action of the quantum 
channel gAB {1® A) gAB , the increase of the von Neu- 
mann entropy AS* for a bipartite n®m state is bounded 
by 

^ S{g%) - ^(pTb) < S{g^Z) - S{g%) + logm, 

(7) 

where S{gA) denotes the entropy of the subsystem A. In 
particular, if the system is separable then AS* < logm. 

Proof.- By the definition the local channel is trace pre- 
serving, hence it does not change the density matrix of 
the first subsystem. Thus p^T* = S^A ^-nd the same holds 
for the corresponding entropies. Then from subadditivity 
of the entropy we have 



S{q7b) < SigT 



) + SigT 



SiQA) + SigT) < 



B 

max 



S{gt) + \ogm. (8) 



We get the first inequality in the Proposition by subtract- 
ing S{g''^g) from both sides of the above inequality. For 
separable states one always has S{g^^^) — Sig'-^) > 0] 
which simplifies (^ to AS' < log m as expected. 

Note that a sequence of quantum channels acting lo- 
cally forms a quantum channel acting locally too. So the 
proposition works also for the dynamics Og ^ and Go, p. 
Moreover, from Eq. (^ we see that the entropy of ini- 
tially pure separable state cannot exceed logm. 
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IV. ENTANGLEMENT VERSUS DEGREE OF 
MIXING 



In this section we study the time evolution of entan- 
glement and compare it with the time evolution of von 
Neumann entropy. To characterise the degree of entan- 
glement we use the entanglement of formation introduced 
by Bennett et al. For any 2(g)2 mixed state this quan- 
tity may be computed analytically as shown by Hill and 
Wootters . In this case the entanglement of formation 
E (or shorter, the entanglement) varies from zero (sepa- 
rable states) to In 2 (maximally entangled states), so in 
the figures we used the rescaled variable E/ln2. 

Our results were obtained by averaging over ensembles 
of random initial states. They were generated according 
to natural measures on: 

(i) 6 dimensional manifold of all pure states for 2 (g) 2 
problem, 

(ii) 3 dimensional manifold of maximally entangled 
pure states, 

(iii) 4 dimensional manifold of separable pure states. 
Numerical experiments have shown that the samples 

of 100 initial states, generated randomly as described in 
the appendix, were sufficient to receive reliable results. 
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FIG. 2. Dynamics of quantum entanglement for system 
Bo,e. Mean entanglement of formation (E) (open symbols) 
and von Neumann entropy (5*) (closed symbols) averaged over 
a sample of 100 random pure states shown as functions of 
discrete time t„. No unitary evolution is present, (a = 0). 
Parameter e, controlling the interaction with environment is 
set to O.Ol(o) or 0.05(D). 



A. Bistochastic channels 

As shown in ||lj,|l5| the mean entanglement of mixed 
states decreases monotonically with increasing degree of 
mixing. Due to interaction with the environment the 
initially pure states become mixed: their von Neumann 
entropy, S{g) = — Tr(gilngi), grows in time. Thus it is 
natural expect a corresponding monotonous decay of the 
mean entanglement. This indeed takes place, as shown 
in Fig. 2 in absence of the unitary dynamics, (a = 0). 
Initial states were taken randomly from the entire space 
of pure states, so in accordance with |T5| , the initial mean 
entanglement is close to (ln2)/2. The parallel processes 
of decay of the entanglement and increase of the entropy 
are accelerated, if the parameter e describing the inter- 
action with environment increases. 



For initially maximally entangled pure states (case (ii)) 
a similar dependence is represented by circles in Fig. 3. 
Here {E{0)) — In 2. The picture changes when unitary 
evolution is involved. The latter leads to oscillations of 
entanglement of formation, reflected in the time evolution 
of entropy. The frequency of oscillations is proportional 
to a. The larger this parameter, the faster the unitary 
evolution U rotates the states g from and into the convex 
set of separable states. In the case of entropy, oscillations 
are only due to changes of the second derivative i.e. en- 
tropy is still monotonically decreasing. This is not the 
case for entanglement E, which can also be seen in Fig. 
4 for several individual initial states (without averaging) . 
For short times the curve for a = (no unitary evolution) 
seems to constitute an envelope for all other curves. 
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FIG. 3. As in Fig. 2 for a sample of 100 maximally en- 
tangled states (S(0) = ln(2)) with e = 0.01; a = O.O(o) and 
a = 0.1(A) for channels described by a) p*-^-* and b) ff^K Ob- 
serve how the influence of the unitary dynamics depends on 
the kind of the channel. 



It is worth to emphasise a significant difference be- 
tween (Fig. 3. a) and 6^^ (Fig. 3.b). In the former 
case the presence of unitary evolution can accelerate the 
process of entropy increase. In the latter, on the con- 
trary, switching on unitary evolution results in slower 
increase of the mean entropy. The iteration of the chan- 
nel Gj, J preserves both the number and the position of 
the nonzero component in p. It is not the case for 9^ ^, 
for which two Pauli matrices generate the entire algebra 
of unitary matrices Ai involved. 
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FIG. 4. Dependence of entanglement of formation on time 
for several randomly chosen maximally entangled pure states. 
The unitary dynamics U = exp{iaH) is governed by the pa- 
rameter a. Here e = 0.01 in and a = 0.1 (narrow lines). 
Reference bold line represents the case of no unitary dynam- 
ics (a = 0), for which the dynamics of entanglement does not 
depend on the initial state. 

Consider now the case (iii), of initially separable states, 
presented in Fig. 5. The presence of the unitary evolu- 
tion may increase the mean entanglement, initially equal 
to zero. However, there is one difference more; for both 
dynamics 9;!^ ^ (Fig. 5. a) and 9^ ^ (Fig- 5.b) presence of 
the unitary dynamics accelerates the process of increase 
of entropy. In absence of the unitary dynamics (a = 0) 
the entropy docs not exceed the value In 2, in accordance 
to our proposition proved in section III. 

Obtained results show that the oscillations of the mean 
entanglement E are anti-corrclatcd with the oscillations 
of the entropy S. It was also checked that if a is kept con- 
stant, but the Hamiltonian is chosen randomly then the 
oscillations of entanglement arc smeared out. It means 
that effects of quantum coherence are destroyed and the 
destruction of entanglement occurs faster. 
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dumped to the ground state. So finally the state of the 
system is a product one consisting of the ground state of 
the affected subsystem and the reduced density matrix 
of the unperturbed subsystem. Thus, after the averaging 
procedure, one gets the averaged von Neumann entropy 
of the subsystem not subjected to action of the channel. 

A random choice of initially pure states of the com- 
posite system induces a certain measure in the space of 
the reduced density matrices |l6j . As proved recently by 
Hall | [l7|] the natural rotationally invariant measure on 
the space of = 4 pure states induces a uniform mea- 
sure in the Bloch ball representing the density matrices 
for N = 2. Denoting the spectrum of reduced matri- 
ces by {1/2 — r, 1/2 -f- r} we may write more formally, 
P{r) = 24r^ for r G [0, 1/2]. The von Neumann entropy, 
averaged over this measure equals 1 /3 [ pT| , in agreement 
with the numerical data presented in Fig. 6. 
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FIG. 5. As in Fig. 3 for a sample of 100 initially separa- 
ble pure states (-B(O) = 0). In absence of unitary dynamics, 
{a — 0), the entanglement equals to zero. 



B. Decaying channel 

Figure 6 presents results obtained for the amplitude 
damping channel (^). In the absence of the unitary evo- 
lution (a = 0) the mean entropy, (S), averaged over the 
entire manifold of pure states (case (i)), does not tend 
monotonically to its maximal value. At t„ 20 the en- 
tropy reaches its maximum and then decreases to its lim- 
iting value about 0.3 (see full circles in Fig. 6b). This is 
due to the fact that for the decaying channel the entropy 
of the system may decrease. 

Numerical data received by averaging over the set of 
maximally mixed states (case (ii), diamonds) and the set 
of separable pure states (case (iii), squares) are shown 
in Fig. 6a. Observe that the steady state limiting val- 
ues of the von Neumann entropy, (S), represents the ini- 
tial average entanglement {E). Indeed, in absence of the 
unitary evolution the perturbed subsystem is eventually 
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FIG. 6. As in Fig. 2 for samples of 100 pure states sub- 
jected to the Kraus channel (4) with p — 0.05 and a = 0.0. 
Initial states drown randomly from the ensembles of a) (□) 
- separable pure states; a) (o) - maximally entangled pure 
states, and b) (o) - ensemble of all pure states. Case (i) with 
unitary evolution, a — 0.1, is denoted by (A) in panel b). 

For non-zero values of a we observe the oscillations 
of the mean entanglement, caused by the unitary evo- 
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lution. It is interesting, however, that the presence of 
unitary evolution allows the final entropy to be maximal 
(see full triangles in Fig. 6b) . It means that the presence 
of the decay channel is completely masked by the unitary 
interaction between both subsystems. 



V. ASYMMETRY OF ENTANGLEMENT DECAY 

We shall consider here dynamics of mixed states hav- 
ing an intriguing property. Namely we choose a quantum 
bipartite system, which violates some entropy inequality 
only with respect to one of both subsystems. Let us recall 
first that the information gain resulting from the mea- 
surement of any of subsystems of a quantum state with 
classical correlations is not greater than the gain obtain 
form measurement performed on the entire system. This 
classical feature is characteristic of quantum separable 
states. They do satisfy the following two inequalities con- 
cerning von Neumann entropy |13,E8|: 



and 



S{gAB) > S^qa), 



S{gAB) > Sign) 



(9) 



(10) 



Shannon entropy, s{x) := —a; In a; — (1 — x) ln(l — x). We 
analysed the time evolution of this quantum system in 
presence of a depolarizing channel 63 ^ given by (||) . In 
the theory of error correcting codes it is one of the most 
popular models of environment induced noise. The evo- 
lution of entanglement for the state p^^^ is represented 
by stars in Fig. 7. In this case the bistochastic channel A 
acts on the 'classical' subsystem B. To investigate a pos- 
sible asymmetry of the entanglement decay we consider 
the state p'^' , for which both subsystems are exchanged. 
More precisely, all elements of both density matrices are 
equal, apart from ^23^ ~ and — P23 ■ The corre- 
sponding dynamics of p^^^ is denoted by crosses in Fig. 
7. In this case the noise interacts with the 'quantum' 
subsystem A. The magnification in the inset reveals the 
asymmetry in the time evolution. Observe that the at- 
tack on the 'classical' part of the system is more harm- 
ful to the entanglement properties of the system. This 
counter intuitive effect, called subsequently anomalous 
entanglement decay (AED), links quantum and classical 
features of the state from information-theoretical point 
of view. 
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where gA and gs denote the reduced density matrices, 
e.g. gA = TrsigAB)- Now we shall focus on the following 
family of states introduced in p^ . They can be written 
gl^Oi^-il + (1 - q)|*2)(*2|, < g < 1, with 



as g 



(1) 



normalised pure state vectors |^i) = a|00) — a^jll) 
and 1*2) = a|10) + \/T~~^|01) with < a < 1. In the 
standard basis, (|00), |01), |10), |11)), the corresponding 
density matrix takes the form 



(l-g)(l-a2) (1 - q)a^J 1 - 

(1 - q)a^J 1 - a2 [1 - q)a^ 





qa ^ 1 — 







q(l - 
(11) 
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Let us take > g > i. Then the first inequality (||) 
is violated, while the second one ( p^ is not. Thus the 
composite system can be called 'quantum' with respect 
to the subsystem A and 'classical' with respect to the 
subsystem B. One may then expect that the bipartite 
system will loose entanglement in different ways, depend- 
ing on whether the environment interacts with classical 
or quantum subsystem. Intuitively one could guess that 
the entanglement should be more robust if the noise af- 
fects the classical subsystem. 

Here we studied the system p^^^ for q = 3/5 and 
=3/4. Then von Neumann entropy of the entire sys- 
tem, 5* = s(2/5) « 0.673 is greater than the entropy of 
the classical subsystem B for which Sb = s(l/4) « 0.562, 
and smaller than the entropy of the quantum subsystem, 
Sa = s(9/20) 0.688, where s stands for the binary 



FIG. 7. Comparison of the dependence of the entanglement 
o: formation for the state g^^^ with = 3/4; q = 3/5 (★) 
aid fj'^-* (x). The bistochastic channel pa with e = 0.01 
ittteracts with the 'classical' subsystem B in the former case, 
and with the 'quantum' subsystem A in the latter case. Solid 
line represents the behaviour of a maximally entangled state 
Pmax- Magnification of the initial dependence provided in the 
inset reveals the asymmetry of the entanglement decay. 



Let us recall that any 2(X)2 system may be described by 
two Bloch vectors, representing locally both subsystems, 
and a correlation matrix T, which represents the projec- 
tion of the composite system onto the family of mixtures 
of maximally entangled states (see [^). A possible ex- 
planation of AED should take into account the fact that 
the local action of environment changes both the Bloch 
vectors, the correlation matrix, as well as their relation- 
ship. A depolarizing channel may affect in a similar way 
both local parameters, but it may distinguish, (in sense 
of the destruction of the entanglement), the correlation 
parameters with respect to the side of the action. 
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It should be noted that, regardless which part is sub- 
jected to the noise, the entanglement of mixed states p''^^ 
and p^^' decreases slower than the entanglement of the 
maximally entangled states (bold line in Fig. 7). This 
is due to the fact that the latter decreases fast for short 
times and slow at longer time scales, for which the ini- 
tially pure state gets mixed. It is thus instructive to 
compare the shape of the bold line starting from In ~ 60 
with the symbols representing the initial decay of entan- 
glement of the states p^^^ . 



VI. AMPLIFYING THE PROCESSES: 
ENTANGLEMENT REVIVALS. 

Consider now the depolarizing dynamics ^ with an 
unitary operation involved, (a ^ 0), affecting either sub- 
systems Aov B. To compare the dynamics of both sym- 
metric mixed states g^^^ and g^"^^ we study their unitary 
interaction governed by the Hamiltonians: H = <^ cr., 



and the reflected one H' = a,, ® a-r 
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Let us consider two cases: 

(a) the noise parameter e is much less than the param- 
eter a characterising the unitary interaction, 

(b) both parameters are of the same order of magni- 
tude. 

Numerical results obtained in the weak noise case (a) 
are presented in Fig. 8. and 9. The revivals of the en- 
tanglement, caused by the unitary interaction, are man- 
ifestly visible, since the strength of the interaction with 
the environment e = 0.002 is much less then the param- 
eter a = 0.06 governing the unitary dynamics. Note 
the characteristic entanglement plateaus, if the analysed 
state travels across the set of the separable states and 
the entanglement attains its minimal value equal to zero. 
The effect of anomalous entanglement decay is clearly 
visible in Fig. 8. a, where the entanglement decays faster 
if the environment interacts with the classical subsystem. 
This contrasts the situation shown in Fig. S.b, for which 
the unitary evolution is due to the reflected Hamiltonian 
H' and the exposure of the 'quantum' subsystem to the 
action of the environment action is more damaging for 
the entanglement. 
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FIG. 8. Unitary dynamics and asymmetry of entanglement 
decay: a) the state g^^^ defined by parameters — 3/4 
and q = 3/5 subjected to the bistochastic cliannel ps with 
e = 0.002 and unitary dynamics H with a = 0.06 (*); the 
symmetric state g^'^^ interacting with the reflected Hamilto- 
nian H' (x). Panel b) shows the data for reflected unitary 
dynamics; the Hamiltonians H' and H are exchanged. 




FIG. 9. As in Fig. 8 for £ = 0.002 and a = -0.06, i. e. the 
process runs bacli in time. Observe that maxima in Fig. 8. a 
correspond to minima in Fig. 9. a and vice versa. 

It is instructive to analyse the same system with the 
unitary evolution reversed in time. Such a case, obtained 
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by a change of the parameter a — > —a, is presented in 
Fig. 9. The general character of the evolution is kept. 
The significant difference is that here the entanglement is 
amplified at the very beginning which may have practical 
consequences if we are interested in short times of the 
process. Note that the figures 8. a. and 9. a reflected along 
the vertical line at tn — (respectively, 8.b and reflected 
9.b) exhibit some kind of symmetry with respect to the 
initial moment. 

What happens if we allow the strength of the coupling 
with the environment to be comparable with the param- 
eter of the unitary interaction? This situation, corre- 
sponding to the case (b), is illustrated in Fig. 10. Here 
some interesting qualitative changes occur. The AED 
effect is present in the case shown in Fig. 10a; at the 
beginning the entanglement disappears faster when the 
'classical' part of the system is affected by the environ- 
ment. Moreover, in this case the entanglement disap- 
pears completely and never revives. If the 'quantum' 
subsystem interacts with the environment, a single en- 
tanglement revival occurs. 
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FIG. 10. As in Fig. 8 for e = 0.01 and a = 0.04. 

In the complementary case, for which g^^'^ interacts 
with the reflected Hamiltonian H' (see 10. b), we observe 
a special kind of competition: for short times the en- 
tanglement is smaller, if the quantum subsystem is per- 



turbed. For longer times, the roles are interchanged, and 
the oscillations of the entanglement are damped faster, if 
the classical subsystem interacts with the environment. 

In general one can see that the pictures correspond- 
ing to the cases (a) and (b) are qualitatively different 
depending on the ratio e/a. This fact may be related 
to the observation concerning the processes of decoher- 
ence. Depending on the relation between two coupling 
parameters the so called pointer basis is determined ei- 
ther by the internal self-Hamiltonian of the system or by 
the Hamiltonian of the interaction with environment . 



VII. DISCUSSION 

We investigated the behaviour of entanglement of bi- 
partite spin-i system subjected to periodic action of the 
environment. The process of destruction of entanglement 
of initially pure states is accompanied by increasing of 
von Neumann entropy. The asymptotic value of the en- 
tropy depends on the form of the interaction with the 
environment. For strongly mixing bistochastic channels, 
(e.g. 8^ and 0'^) the entropy achieves the maximal value 
In 4. If the decaying channel is involved, the entropy gets 
its maximum and then it monotonically decays to the 
asymptotic value, which reveals the initial entanglement 
of the system. 

If the internal unitary evolution entangling the system 
is present, the decay of the entropy due to the decaying 
channel can be replaced by the process of mixing the state 
more and more. The general feature of the time evolution 
is that the entanglement decreases as the system becomes 
more mixed. This corresponds to the results recently 
obtained in ||lj,|l^, where it was shown that the mean 
entanglement of quantum states, averaged over a sample 
of mixed states with the same von Neumann entropy, 
decreases with the degree of mixing. The presence of 
the internal unitary evolution leads to the revivals of the 
entanglement and to suppression (or acceleration) of the 
entanglement decay. 

Perhaps the most intriguing is the character of asym- 
metry of the time evolution of the entanglement. For 
some initial mixed states consisting of two non-equivalent 
subsystems, the entanglement decays faster, if the 
environment interacts with the 'classical' subsystem, 
which satisfies the entropic inequality. Many years ago 
Schrodinger considered entanglement of pure state as a 
property of having both subsystems less informative for 
the observer, than the composite system. Mixed states 
( [III ) considered here exhibit this property only with re- 
spect to one subsystem . Our results show that the ac- 
tion of environment to the 'classical' subsystem is some- 
times more harmful to the entanglement. In this case one 
can thus say that the quantum entanglement runs away 
faster through the classical door. 

In the context of the above discussion some general 
questions emerge. Consider a quantum entangled state g 
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with, say, S(gA) > S{gB), not necessarily violating the 
inequality (^j . Under which conditions the entanglement 
is less robust to action of the environment on this sub- 
system, for which the entropy of the reduced operator is 
smaller? How is it related to a possible violation of the 
von Neumann entropy inequality by subsystem A7 What 
happens if instead of the inequalities (^,^|) one applies 
the generalised a-entropies inequalities ( satis- 
fied for classical systems? All these questions seem to be 
important for deeper understanding of the dynamics of 
quantum entanglement. 

It would be also interesting to analyse the role of en- 
tropic asymmetric states like (^Tj) in context of quantum 
communication. In fact these states have only one co- 
herent information positive |^ , p3| (see also ||2^). For 
the corresponding quantum channels this might imply an 
asymmetry in the transfer of quantum information with 
respect to its direction, {A ^ B oi B ^ A). 

Finally, the obtained results show that even the sim- 
plest bipartite systems may exhibit non-trivial properties 
form the point of view of the information theory. In this 
context it would be important to investigate further the 
dynamics of mixed entanglement, in particular, by tak- 
i ng^ into account the phenomenon of bound entanglement 
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APPENDIX A: RANDOM PURE STATES 



I*) = (cosi93,sin^l3COSz?2e''^', 
sin i?3 sin i?2 cos -die^f^ , sin -ds sin '02 sin z9ie"^i ) , (Al) 

where -dk € [0, 7r/2], and ipk € [0, 27r) for fc = 1, 2, 3. 

A uniform distribution over almost all of CP'^ is ob- 
tained by choosing the uniform distribution of the 'az- 
imuthal' angles; P{fk) — l/27r. In the analogy to the 
volume element on the sphere the 'polar' angles should 
be taken in a nonuniform way, with the probability den- 
sity ^ 

P{dk) = ksm{2^k){sm^kf''~^ (A2) 

for dk G [0,7r/2], k = 1,2,3. In practice it is convenient 
to use auxiliary independent random variables dis- 
tributed uniformly in [0, 1] and to set dk = arcsin(^^^^'') . 
Above formula with fc = l,2,...,A^— 1 allows one to get 
a natural distribution on CP^^^ [G7|. 



2. Random separable pure states 

Any 2® 2 pure separable state may be written as \^s) — 
I'ipi) ^ IV'2), where \tpi) and 11/12) are N = 2, one-particle 
pure states. The 4 dimensional manifold of separable 
states has thus a simple structure of a Cartesian product 
CP^ X CP^. a uniform measure on this manifold is ob- 
tained be taking both states ji/'i) distributed uniformly 
(and independently) at the Bloch sphere, CP^ ^ S^. 

Working in the standard basis, 

I*.) =C/i®t/2|(l,0,0,0)), (A3) 

where Ui and U2 denote two independent random uni- 
tary matrices distributed uniformly on SU{2). This 
parametrisation describes the entire 4D manifold of the 
separable pure states. 



In this appendix we present algorithms allowing one to 
generate random quantum states distributed uniformly 
at the entire space of pure states, the manifold of sepa- 
rable pure states and the space of maximally entangled 
pure states. We concentrate here on the simplest 2 (g) 2 
problem, but the algorithms below can be easily gener- 
alised for higher dimensions. 



1. Generation of random pure states 

The set of pure states of a 4 dimensional Hilbert space 
forms a complex projective space CP'^, on which a nat- 
ural, unitarily invariant measure exist. To generate ran- 
dom pure states according to such a measure on this 6 
dimensional space we take a vector of a random uni- 
tary matrix distributed according to the Haar measure 
on C/(4). The Hurwitz parametrization pffl gives 



3. Random maximally entangled states 

In an analogous way we may represent the maximally 
entangled states as 



I®i7i|(0,l,l,0)/\/2). 



(A4) 



It is easy to see that for this states the reduced den- 
sity matrix is proportional to identity matrix, and the 
entropy of entanglement achieves its maximum In 2. The 
states obtained by a symmetric operations Ui®l. are also 
maximally entangled. Using the standard representation 
oiUi we parametrise maximally entangled states by pq] 



l*el) 



1 
71 



cosiJe^'^i 
-sini?e-*'^2 



(A5) 



9 



The angles are di strib uted uniformly in [0, 27r), 
whereas according to (Al) = sin(2z9) for -d e [25] 
[0, 7r/2]. Note that the standard element of the vol- 
ume on the two sphere dS = sinOdOdip is written in a 
rescaled variable 6 = 2'd. Given maximally entangled [26] 
state corresponds to a single unitary matrix Ui pertain- 
ing to SU{2), but the 3-D manifold of the maximally [^'''l 
entangled states has the topology of the real projective 
spaq o, W.P^ r-,u{3)/Uil) M. [281 
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FIG. 1. Discrete model of periodic dynamics (cf. Fig. 
8.1, 8.2 of Ref. 0]). Interaction with the environment A 
transforms the state Qn into q'„ and then the unitary trans- 
formation U maps it into Qn+i- 
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